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Abstract 

First principles computations of Raman intensities were performed for perovskite-family com- 
pound CaAL^NbirtOa (CAN). This compound features 1:1 (NaCl-type) ordering of Al and Nb 
superimposed onto the b~b~c+ octahedral tilting. Raman tensor for CAN was computed using 
the package for first-principles computations ABINIT (URL http://www.abinit.org). Computa- 
tions performed for both untilted cubic (Fm3m) and tilted monoclinic (P2i/n) CAN structures 
showed that the strongest Raman lines are associated with the ordering of Al and Nb. The com- 
puted spectrum agreed qualitatively with the experimental data measured on powder (CAN is 
available in polycrystalline form only). The effect of cation disorder on the Raman intensities was 
considered using phenomenological theory of light scattering in the vicinity of a phase transition. 
We suggest that, for certain modes, the corresponding Raman intensities depend primarily on the 
average long range order while, for other modes, the intensities are determined by fluctuations of 
the order parameter. 
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I. INTRODUCTION 



Electronic properties of double perovskites AB / 1 / 2 B^ 2 03 are strongly influenced by the 
degree of B-cation ordering.i^ Raman spectroscopy is an attractive technique for probing B- 
site ordering in perovskites.^ Raman spectra proved to be sensitive even to small nanoscale 
ordered regions. For example, Raman spectra of largely disordered PbSci/2Tai/203 still 
feature well detectable lines attributed to the existence of nanoscale ordered regions. 5 

The quantitative theoretical description for Raman spectra of disordered or partially 
ordered double perovskites has not been developed. Furthemore, no first-principles compu- 
tations of Raman intensities in completely ordered double perovskites have been reported. 
This lack of theoretical analyses of Raman intensities impedes development of Raman spec- 
troscopy into a quantitative tool for probing B-cation ordering in perovskites. 

Raman effect refers to a light scattering by normal vibrations having appropriate sym- 
metry. Landau and Placzek^ have shown that Raman intensities can be expressed over the 
phonon-mediated change of high-frequency (electronic) dielectric permittivity s{ujl) where 
ujl is the incident light frequency of a laser. The theory commonly exploits this relation to 
evaluate the effect of a specific vibration on Raman spectra. This approach requires complex 
computations of the high-frequency (electronic) dielectric permittivity tensor and its weak 
variations caused by ionic vibrations. First-principles computations of Raman spectra have 
been limited to several semiconductors (Si, Ge, SiC>2 etc.) and rare-earth transition-metal 
borocarbides . 7l8l9llQl11 ) 

The present study reports the first attempt to compute Raman intensities for double 
perovskites from first principles. The paper is structurally as follows. Section 2 describes 
the general theory for the computation of Raman intensities. Section 3 considers ordered 
double perovskites AB'^B^Oa and includes three parts: the general theory of Raman 
spectra for cubic double perovskites (Section 3.1), results of the computations for the untilted 
cubic CaAli/ 2 Nbi/203 (CAN) (Section 3.2), and results of the computations for the tilted 
monoclinic CAN (Section 3.3). The effect of B-cation disorder on the Raman spectra of 
double perovskites will be described in Section 4. Section 5 considers the effect of phonon 
dispersion on the Raman line profile. Finally, Section 6 summarizes the results obtained. 
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II. THEORY 



Born an Huang 12 derived an expression, which reduces the computation of integral Raman 
intensities to the evaluation of the imaginary part of linear (Raman) susceptibility (we 
consider the part of this intensity connected with the mode i, polarization directed along a, 
and the field along (5): 



I iaf s = K,R ia p J [n (u) + 1] Imxia/3 {oj)du (1) 

— oo 

where k = 2h (ujl — Ui) 4 / (vrc 4 ), lul is the laser frequency, Ri a p is Raman tensor, n(uj) + 1 = 
[1 — exp (— hu i /k B T)}~ 1 is Bose occupation number, Ui is the frequency of mode i, k B is 
Boltzmann constant, h is Planck constant, T is temperature, and Xia/3 is the contribution of 
mode i to dielectric susceptibility. Kramers-Kronig transformation yields^: 

hap = 2H ( ? " ^ [n (co,) + 1] {R iaP f (2) 

Thus, the relative intensity of the Raman active modes is determined by the fundamental 
characteristics of the lattice dynamics: phonon frequency, mass, and Raman tensor. 

The Raman tensor is a derivative of the electronic part of the dielectric permittivity 
tensor e a p at a laser frequency with respect to the z-th normal mode displacement: 

Riaf, ~ (3) 

The first-order Raman scattering can be described using the expansion of polarization 



Pa — PaO + £0 X] P-iapEpQi — 
i 

= P a0 + e ' P4ap, ni Epe ini m~ 1,2 Qi (4) 



where e m7 is the eigenvector of a dynamical matrix, m n is the mass of the atom n, and Qi 
is the normal mode displacement magnitude. From Eq. 



n-y 
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Eq. (J3J) can be used to estimate the contribution of different ions in mode i into the Raman 
tensor. 

The integral Raman intensity can be expressed as: 

2h (u L - u.f 2 
Uap = 7 n (Wi) + 1 (r ia/3 ) (6) 

with 

Vm / 

where mo = la.m.u. 

Expression (J2J) is applicable to single crystals. For powders, the integral intensity has 
be averaged over the possible directions of crystallites which can be accomplished using the 
rotation invariants 



2 (j'ixx ^iyy 7'izz) 



Vizy ^iyz) 



G? 2 [(r ixy + r iyx f + (r txz + r lzx f + 

(jizy Tiyz) ] H l~g [(''in ^iyy) 
~\~ (j'ixx fizz) (j"iyy fizz) ] 



(8) 



For the 90° and backscattering geometry, a (polarized) laser beam, and a scattered radiation 
analyzer in use, the (reduced) intensity for the powder can be expressed a o 6 ' 14 ' 15 : 



T powder / \4 l4-n(u>i) 

jpowder ^ _ ^4 

powder jpowder , ypowder 



10Gf } + 4GP 
5GP + 3GP 



(9) 



jpowaer jpowaer , j 

hot — 1 w + 



powder jpowder / jpowder 

Pi — i _L 



where subscripts || and _L correspond to the polarized and depolarized light, respectively. 
Finally, one can plot the Raman spectrum by using the expressions: 
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jpowder / \ v-> rpowder Fj 



oowder J\ 



jpowder / \ ^ T powder r, / ln \ 



rpowder / \ v— > rpowder V 4 

where Tj is the damping constant for the i-th mode. The computation of this constant should 
be done separately. Here we concentrate on the computation of integral Raman intensities 
only. 



III. RAMAN SPECTRA OF ORDERED AB' 1/2 B^ /2 3 PEROVSKITES 

A. Results of first-principles computations for cubic CaAlo.5Nbo.5O3 

CaAlo.5Nbo.5O3 (CAN) crystallizes with a perovskite-like structure and exhibits 1:1 or- 
dering of Al and Nb superimposed onto the b~b~c + tilting of oxygen octahedral This 
combination of cation ordering and octahedral tilting yields monoclinic P2 1 /nsymmetry 
with 20 atoms per unit cell.— In the absence of tilting, the structure would adopt a cubic 
Fm3m symmetry yielding A\ g + E g + 2i<2 9 Raman-active modes (see Fig. and Ap- 
pendix). First-principles computations of vibration frequencies, dynamical charges, and 
infra-red spectra for the cubic and tilted monoclinic structures of CAN have been reported.— 
The present study focuses on the computation of Raman intensities with the help of a package 
for ab-initio computations, ABINIT, which is a result of joint project between the Universite 
Catholique de Louvain, Corning Incorporated, and other contributors. 18 Ingredients needed 
for the calculation of responses to atomic displacements and homogeneous electric fields have 
been already described in the literature.— 1 ^ We do not compute damping constants in the 
present study, only integral intensities. A damping constant Tj = 10 cm -1 was used for all 
modes. 

First, we performed the first-principles computations for the 10-ion unit cell of CAN 
having symmetry Fm3m (no tilting). We used 8x8x8 Monkhorst-pack grid, ECUT=45 
Ha and tolvrs = tolwfr = 10~ 18 . We computed high frequency (electronic) dielectric tensor 
for the 5 different normal mode displacements. Then, we interpolated the results with a 
polynomial and extracted the linear dependence of the high-frequency (electronic) dielectric 
permittivity on the normal mode displacement. The results of the computation of the Raman 
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tensors are summarized in Table 1. The A\ g mode has the largest Raman tensor. Table 2 

presents the coefficients needed for the computation of the average Raman tensor for the 
F 2g modes. Results of the computation of Raman intensities are presented in Table 3 and 
Fig.H 

B. Raman spectrum of tilted CAN structure (P2i/n) 

For the tilted CAN (P2i/n), we used the eigenfrequencies and eigenvectors which are 
similar to those obtained by Cockayne. 1,6 The lattice constants in these computations were 
taken from the experiment, but the distortion and the equilibrium atomic coordinates were 
determined selfconsistently under the constraint of a fixed volume. The obtained structure 
contains frozen tilting of the b~b~c + type in the Glazer notations, and the equilibrium 
structure is monoclinic P2 x /n with the unit cell 4 times larger than the primitive perovskite 
cell. By using ABINIT, we computed high-frequency (electronic) dielectric tensor, £ a p(^L), 
for each of the displacements Q. In all cases, we used the same high precision and 8x8x8 
Monkhorst-pack grid as in the computation for the cubic CAN. To facilitate the comparison 
with experiment, we used the damping constants of 8 cm -1 and 2 cm -1 at frequencies above 
and below 400 cm -1 , respectively. 

Fig. presents the experimental (a) and the calculated total (b), polarized (c), and 
depolarized (d) Raman spectra for a ceramic CAN (the integral intensities for all lines are 
shown in Table 5). 

IV. DISORDERED AB' 1/2 B^ /2 3 PEROVSKITES 

Raman spectra of double perovskites reportedly are sensitive to the degree of cation order, 
The 1:1 order parameter for the Pm3m — > Fm3m transition is a scalar and transforms 
according to the Rf representation of the parent Pm3m space group. A useful insight into 
the dependence of Raman intensities on r\ can be obtained using phenomenological theories 
of light scattering in the vicinity of a phase transition . 13 > 21 We start with the expression 
(J7j)and expand the Raman tensor into a series with respect to the order parameter: 

r ia p{r) = r iaPtl r](r) + ... (11) 
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The change of dielectric permittivity e(r) is proportional to the product of 77 (r) and 
displacement Q(r). The Fourier transform of this product can be expressed as: 

A^~5>k +q Q-k (12) 

k 

Finally, the intensity of the light scattering is proportional to the statistical average of the 
square of |A£ g |: 

1 ~ (^k+qQ-k^-k'-qQk') (13) 

kk' 

According to this expression, phonons in the entire Brillouin zone are excited if \r)\ < 1. 

Recent experimental studies of ^-dependence of Raman intensities in CAN have shown 
that the intensities corresponding to both F 2g - and E^-related modes vary approximately 
as square of the order parameter.— In contrast, the intensity of the A\ g line exhibited weak 
dependence on 77. These observation can be understood by considering two particular cases 
of the general expression (JT3J). 

First, we consider the case when Q does not correlate with 77: 

~ E ( W?-k- q > (QkQ-k> (14) 

k 

Now, we present 77 (r) as the sum of the average order parameter, t] e and its fluctuation 
Ar/(r): 



~ vl (QnQ-n) + 
+ £ (Ar/ k+q A77_ k _ q ) (g k g_ k ) (15) 

k 

The second term is finite only in the vicinity of the phase transition temperature, and it is 
finite even above T c (the k = contribution behaves^ as |r| -1 / 2 ). The first term behaves 
as: 

T —T 

V 2 e -r= I ^-A (16) 

1- a 

where T c is the order - disorder phase transition temperature, T a the annealing temperature 
(the equation can be applied to the samples quenched from T a ). This type of the dependence 
was experimentally observed for the F 2g and E g related modes in CAN.— 
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Now we consider the case when Q and rj are strongly correlated and the displacements 
are modulated by the order parameter: 

Q(v)~ V (r) (17) 

In particular, the displacements change sign at the antiphase domain boundaries following 
the change in the sign of the order parameter, so that changing the sign of the order 
parameter has no effect on the Raman intensity Even at the antiphase domain wall, where 
1] varies linearly from -1 to 1, one gets a finite contribution to the intensity. Relation (JT7jl 
is realized only for those modes which transform similarly to the order parameter. In the 
symmetry group, Fm3m, only the Ai g mode transforms similarly to Rl + mode in Pm3m. 
In this case, fllHj) takes the form: 

# = Ctf e (A^_ q ) + Drf e + 

+Ei £ {Arj cl+K Ari_ K Ari_ (l+v Ari_ v )+ (18) 
+Fii] e £ (A?? q+K A?7_ K A?7_ q )... 

where^i 

(lA^I 2 ) ~ r- 1 (19) 

The first term is temperature independent below T c and should vanish above T c . The second 
term is negligible in the vicinity of a phase transition (it varies as r 2 ), and the last two terms 
are spread over phonon bands (they are finite above T c ). The third term depends 2 ^ on r as 
|r| -1 / 2 ). Experimental observations^ 2 - suggest that, the main contribution to light scattering 
for the A lg mode comes from the fist term in (39) and is temperature independent. 

Thus, we conclude that the F 2g and E g modes obey dependence (|T7)|) whereas A\ g mode 
shows mostly constant intensity below T a = T c . This difference stems from the fact that the 
Aig mode transforms in the same way as the chemical order parameter while the F 2g and 
Eg modes transform differently. 

The F 2g and E g modes interact with strains, which can be constructed from the com- 
ponents of the gradient of the order parameter. This interaction mixes the polarized and 
depolarized spectra. 

In the distorted (tilted) structure of CAN, there are Raman lines originating presumably 
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from tilting but not from the chemical ordering. A ferroellastic order parameter should 
supplement rj, in this case. 13 

V. PHONON BAND WIDTH 

Disorder can cause excitation of the g-points in the entire Brillouin zone . 23 ' 24 ' 25 ' 26 The 
contribution of this effect to the width of the Raman line is limited by the width of the 
corresponding phonon band. We have performed first-principles computation of phonon 
frequencies for a 80-ion cubic supercell of CAN (Viena DFT package VAS P 27 ' 28 was used) 
and compared the results with those obtained for the tilted 20-ion supercell (Fig. EJ). 
According to our results, the upper branch of the phonon bands (which includes the point 
corresponding to the A\ a mode) has the width of about 70 cm -1 , which is comparable to 
the width of the A\ g Raman line detected experimentally (about 40 cm -1 ). Such narrow 
width increased only slightly with temperatures approaching T c , which can be attributed 
to relatively narrow width of the corresponding phonon band. This relatively narrow width 
of the high frequency phonon bands in CAN can be caused by the band folding due to a 
doubling of the unit cell upon ordering. 

VI. SUMMARY 

We have performed ab initio computations of Raman intensities for the AB^B'/^Os 
perovskites. The computation of the Raman intensities in cubic CAN revealed that all 
Raman-active lines, A lff +E 9 +2F2 g , exhibit large relative intensities in single crystals. The 
low-frequency part of the Raman spectrum computed for the untilted cubic CAN structure 
differs significantly from the corresponding part of the Raman spectrum for the tilted mono- 
clinic CAN. The intense lines in the cubic structure at 101 cm -1 shift to higher frequencies, 
split (188 cm -1 , 242 cm -1 , and 211 cm" 1 ),— and become relatively weak. This is a result 
of the change in frequency, occupation number, and Raman tensor for this line. The high- 
frequency F2 g line in the cubic structure at 445 cm -1 splits in the low-symmetry structure 
into three lines at 446 cm -1 , 452 cm -1 and 457 cm" 1 .— The relative intensities of these lines 
decrease in the distorted structure though to a lesser extent than for the low-frequency F2 9 
line. E g line at 660 cm"" 1 splits into two lines at 553 cm -1 and 554 cm" 1 ,— and the integral 
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intensity for this line decreases. Ai g line at 910 cm -1 in the cubic structure exhibits signif- 
icant downshift to 797 cm -1 in the distorted CAN.— Experimentally, this line is observed 
at higher frequency of 844 cm -1 . This discrepancy between the experimental and computed 
frequencies can be due to the overestimated degree of tilting in the DFT computations. The 
computed Raman intensities agree qualitatively with those observed experimentally. Yet, 
the experiment yields relative integral intensity of the K\ g line larger than the computed 
value (mostly due to the larger width), and the computed E g intensity is smaller in the 
computation in comparison with experiment. This can be related to the contribution from 
the scattering on the powder surfaces. Second order Raman processes can contribute to the 
detected line intensity as well. 

On the basis of the present first-principles computation we conclude that: (i) In agree- 
ment with earlier speculations,— high-frequency Raman lines appear in AB'^B'^Os dou- 
ble perovskites mainly due the ordering of B' and B" ions; (ii) The modulation of the 
ionic vibrations by the order parameter in CAN can explain the weak dependence of the 
A\ g line intensity in CAN on the annealing temperature when approaching T c ; (iii) Our 
first-principles computation has shown that the width of the high-frequency phonon band 
is about 70 cm -1 that can explain the weak dependence of the high-frequency Raman lines 
asymmetry in CAN on the degree of disorder. 
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APPENDIX A: GENERAL EXPRESSIONS FOR CUBIC Fm3m STRUCTURE 

Cubic Fm3m ordered perovskite AB'^B'^^Oa exhibits Raman active modes Ai g + E g + 
2F2 g . Ai g is a fully symmetric breathing vibration of oxygen octahedra (Fig. Q). The Raman 
tensor for this mode is diagonal, 
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is 



a 



a 
\0 a J 



(Al) 



with the diagonal element 



a = r A 



(6m Q ) 



'1/2 



?" z2,Ol£ — r zz,02z + J" zz,03x~ 



T zz,OAx f zz,05y ^zz,06y) 

here the oxygen atoms in an octahedron are numbered from 1 to 6 (Fig. 
Symmetry of the lattice requires that: 



(A2) 



T"zz,02z T"zz,01z T"xx,03x Txx,04x 

Tyy,05y ^"yy,06y 
T"zz,03x ^zz,04:X ^zz,05y 'zz, 06y 

^xx,01z Txx,02z ^xx,05y 'xx,Q6y 

— fyy,Olz — —fyy,02z = fyy,03x = —fyy,OAx 



(A3) 
(A4) 



Hence, (IA2I) reduces to: 



m 



i r zz,oiz + 2r 22 



03x) 



(A5) 



where mo is the mass of an oxygen atom. Thus, the Raman tensor for this mode is deter- 
mined by the two quantities describing the longitudinal and transverse changes of polariza- 
tion with respect to atomic displacements. 

Eg mode corresponds to another type of breathing of oxygen octahedra (Fig. Q). This 
mode is doubly degenerate. The corresponding Raman tensors can be represented as: 



( -b \ 



r E 3 m 



-6 
2b 



( c (A 
-c 




(A6) 



where 
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-(2r 



zz,01z 



- 2r 



4 v / 3mo 

T"zz,03x ^zz,OAx T"zz,05y ^zz,06y) 
1 



zz,Q2z~ 



(A7) 



r Eg (2)z 



\fzz,02,x — 



— r zz,OAx ~ r zz,05y + r zz,06y) 

From the symmetry relations (jA3|) : 



b 

c 



i 



1 



(jzZ, Ol2 ~~ r 2^,03xJ 

0" zz.Olz — r zz,03x) 



(A8) 



Therefore, the Raman tensor for the E g modes is controlled by the same two constants as 
the Raman tensor for the A\ g mode, and 



y/3b. 



(A9) 



F2 9 modes are triply degenerate and contain contributions from the oxygen and A-cations. 
Oxygen-related F2 9 o modes are presented in Fig. d Their Raman tensors are nondiagonal: 



where 



r F 29 o|2J 



Xo0 



d 
d 
\0 J 



d Q \ 




V^o J 

/oo o ^ 

d 

\o d o y 



(A10) 



V xy,03x ^xy,OAx 



(All) 



^ 2^7^ v ' xy ' 

Because of the symmetry, only two ions out of six (those residing on axis x) contribute to 
this Raman tensor. The symmetry relations give: 
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xy,03x 



i r 'xy,Oix T"xz,03x 



T"xz,04:X 
T"yx,05y 
~ r zx,Q2z 



Tyz,Qhy 
~T"yx,06y 
^"zy,01z 



^yz,OQy 
f zx,01z 
— r zy,02z 



(A12) 



Thus, only one constant determines the Raman tensors for the F 2g o modes: 



o 



'mo 



Txy,03x 



(A13) 



The A-cation related (F 2g A) modes (Fig. ^) have the same Raman tensor as the oxygen 



related F 2g modes: 



( d A \ 
rf A o 

\o o o y 

o o 4^ 
ooo 

\d A J 

^ o ^ 

d A 

V o d A o y 



where 



y/2m A 



^xy,Alz ^xy,A2z) 



with the symmetry relations: 



(A14) 



(A15) 



'"asj/jAla ^'xy,A2z ?":rz,aly 

Txz,A2y ^yz,Alx 
= r yz, A2x = r yx,Alz = r yx,A2z 
Tzx,Aly T zx,A2y 



From these relations: 



(A16) 
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rf A = r F2gA $ >xy = y ^r xy , Alz (A17) 

The eigenvectors of the F2 5 modes are linear combinations of the O- and A-related normalized 
vectors: 



e F2g = c e F2g0 + c A e F2gA (A18) 

Their Raman tensors can be represented as linear combinations of the O- and A-related 
Raman tensors 

^F 2ff a/3 = C r F2g0a p + C A r F2g Aa/3 ( A19 ) 

The polarized I^ yst and depolarized ij C _[ ys * Raman spectra of the cubic AB^/ 2 B^ 2 03 P er_ 
ovskites are determined by the diagonal and nondiagonal elements of the Raman tensor, 
respectively: 

rcryst r \4 r 



ir^K-^r^f [l + n( Wi )] 



jcryst 



(u L - LOif &*.[l + n fa)] (A20) 



jcryst jcryst . jcryst 

From these quantities, the corresponding Raman spectra can be computed as: 



J yst (oj) 


- . v T cryst Ti 




-cryst ^ 


T cryst Ti 


(A21) 


-cryst ^ 


v - ^ jcryst Ti 





The Raman spectra for powders can be found from expressions (jl(Jj) . 
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TABLE I: Relative values of Raman tensor components obtained from the first-principles compu- 
tations of cubic CAN. 



notation 


value 


a 


0.03084 


b 


0.00717 


do 


0.01572 


dCa 


0.00375 



TABLE II: Eigenvectors for Fig modes. 



coefficients 


F 29 I 


F 23 H 


co 

CCa 


0.9794 
-0.2019 


0.2019 
0.9794 



TABLE III: Relative intensities of Raman lines obtained in the first-principles computation for 



single 


crystal cubic CAN (T 


= 300K). 










mode 


v [cm -1 ] 


jcryst 


Geometry 


jpowder 
1 t 


jpowder 

h 


jpowder 


A lg 


910 


8.7388 




8.7388 


8.7388 


0.0000 


E 9 


660 


2.8332 




1.6527 


0.9444 


0.7083 


F 23 I 


445 


4.9791 


_L 


6.9708 


3.9833 


2.9875 


F 2ff H 


101 


11.8411 


_L 


16.5776 


9.4729 


7.1047 
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FIG. 2: Computed Raman spectrum of cubic CAN: I? yst (a), I^ yst (b), I c [ yst (c), lf owder 

jpowder / \ jpowder /r\ 
1 \\ \ e )i W- 
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FIG. 3: Computed Raman intensities for ceramic CAN (P2i/n): experimen t 17 ! 22 (a), theory: 

T powder \ T powder / \ T powder / n \ 

n (b), i\\ (c), r ± (d). 
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FIG. 4: Results of first principles computations of phonon frequencies in different supercells of 
CAN: the data for the 20-ion supercell^ are shown for comparison with our computation for 
80-ion supercell of CAN. 



TABLE IV: The relative integral intensities obtained in the computation of CAN (P2\/n). The 
stars show the modes, which are closest to those obtained under the Fm3m constraint. 



v [cm 


jpowder 


jpowder 


jpowder 


797* 


1 


2034 


0.0013 


1.2048 


554* 


o 


0114 


0086 

V / • V / V / V.- 1 V ' 


0200 


553* 


n 


0149 


0.0100 


0.0249 


457* 


n 


1 927 


0.1442 


3368 


452* 


o 


1698 


0.1273 


0.2971 


446* 


o 


1443 


0.1082 


0.2526 


242* 


o 


0306 


0.0228 


0.0534 


211* 


n 


01 69 


0.0127 


0296 


188* 


o 


0324 


0.0241 


0565 


770 


o 


0004 


0.0003 


0.0008 


351 


o 


0015 


0.0011 


0026 


732 


o 


0027 


0.0020 


0.0047 


482 


o 


0179 


0.0134 


0.0312 


289 


o 


0000 


0.0000 


0.0000 


331 


f) 


0168 


0.0113 


0.0281 


237 


o 


0072 


0.0054 


0.0126 


359 


o 


0113 


0.0084 


0.0197 


564 


n 


0030 


0.0022 


0052 


185 





0026 


0.0019 


0.0045 


175 





0161 


0.0121 


0.0282 


136 





0109 


0.0077 


0.0186 


182 





0018 


0.0014 


0.0032 


162 





0063 


0.0042 


0.0105 


339 





0079 


0.0025 


0.0104 
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